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Binary black holes can form efficiently in dense young stellar clusters, such as the progenitors of 
globular clusters, via a combination of gravitational segregation and cluster evaporation. We use 
simple analytic arguments supported by detailed iV-body simulations to determine how frequently 
black holes born in a single stellar cluster should form binaries, be ejected from the cluster, and 
merge through the emission of gravitational radiation. We then convolve this "transfer function" 
relating cluster formation to black hole mergers with (i) the distribution of observed cluster masses 
and (ii) the star formation history of the universe, assuming that a significant fraction g c \ of star 
formation occurs in clusters and that a significant fraction g e va P of clusters undergo this segregation 
and evaporation process. We predict future ground-based gravitational wave (GW) detectors could 
observe ~ 500((? c i/0.5)(<;ovap/0.1) double black hole mergers per year, and the presently operating 
LIGO interferometer would have a chance (50%) at detecting a merger during its first full year 
of science data. More realistically, advanced LIGO and similar next-generation gravitational wave 
observatories provide unique opportunities to constrain otherwise inaccessible properties of clusters 
formed in the early universe. 
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Given our understanding of how isolated binary stars 
evolve, noninteracting stellar systems should produce rel- 
atively few double black hole (BH-BH) binaries tight 
enough to merge through the emission of GW within 
the age of the universe [l], 0, HI- Portegies Zwart and 
McMillan 0] demonstrated that interactions between 
black holes (BHs) in dense cluster environments could 
produce merging BH-BH binaries much more efficiently 
than through the evolution of isolated binaries. As a re- 
sult, the local binary black hole merger rate - the net 
rate both from isolated evolution of noninteracting stars 
and from dense clusters - can depend sensitively on the 
formation and evolution of young clusters through the 
entire history of the universe. 

An increasing number of galactic [B|, 0] and extragalac- 
tic @,[|| observations suggest at least 20% of stars form in 
dense, interacting stellar clusters. Over time, each cluster 
dissipates, both because hot young stars and supernovae 
(SN) heat and eject a significant fraction of the residual 
gas that gravitationally binds the cluster ("infant mor- 
tality"), and because the host galax y's tidal field strips 
off stars as the cluster orbits it 0, Il0l | . Thus any set 
of coeval clusters decrease in number and size, spewing 
stars into their hosts [ll|, [l2| , with only a few of the most 
initially dense and orbitally-fortunate clusters surviving 
to the present. 

Unfortunately, electromagnetic observations of large 
young clusters in other galaxies cannot resolve their inter- 
nal structure. These observations therefore only weakly 
constrain the fraction of young clusters that survive their 
first few Myr, during which the most massive young stars 
evolve, supernovae, and give birth to black holes. Simi- 



larly, they cannot determine the efficacy of runaway stel- 
lar collisions, in which massive stars in particularly dense 
clusters gravitationally segregate and collide to form very 
massive stellar progenitors to large 100 M Q ] single 
or binary black holes 0, 0, EE E3]- How- 
ever, clusters give birth to many binary black holes, 
whose GW merger signals provide unambiguous infor- 
mation about the processes which produced them. For 
example, Fregeau et al. [l3| demonstrated that runaway 
stellar collisions produce very massive [~ 100 M@\ dou- 
ble BH binaries often enough to be easily and unam- 
biguously seen with advanced ground-based interferom- 
eters (i.e., advanced LIGO/VIRGO). The detection rate 
of such massive BH mergers would therefore constrain 
the fraction of young clusters which undergo collisional 
runaway. Similarly, in this Letter we argue that a high 
GW detection rate of ~ 15 — 20 M Q black holes will un- 
ambiguously measure how often young clusters dynami- 
cally evaporate their most massive compact components. 
Specifically, we use the results of a recent set of numerical 
simulations of BHs in clusters to present a simple, 
analytic approximation relating the BH-BH merger rate 
to properties of young clusters. Existing observations 
only weakly constrain the first few Myr in the lives of 
the most massive and ancient clusters; advanced ground- 
based gravitational wave detectors therefore can explore 
new facets about the evolution of these first massive clus- 
ters. 

Within each cluster, stars form accor ding to an ap- 
proximately power-law mass distribution [19| , with only 
a small number as 3 x 10 _3 (M c ;/Mq) of stars more 
massive than 20M Q , roughly the mass needed to guar- 
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antee a black hole forms in a SN . These most mas- 
sive stars should rapidly evolve, undergo SN , and form 
black holes on a timescale much shorter than the relax- 
ation time of the cluster. The most massive clusters with 
M c i > M cr it = 3 x 10 4 should contain more than 100 such 
stars, each of which produces a roughly rribh, ~ 10 — 25Mq 
hole [20j. Being the most massive objects in the cluster, 
systems (single or binary) containin g B Hs will rapidly 
mass segregate to the cluster core [H, HH • In the process 
exchange interactions will quickly break up any remain- 
ing star-BH binaries Since the numerous black 
holes significantly outweigh the average constituents of 
the cluster ((m*) « 0.5M Q ), the black holes can decouple 
from the stars in a process known as the "Spitzer instabil- 
ity" [23], [3] ■ The more massive BHs interact and evolve 
on a more rapid timescale ~ t c i (to*) / {rribh) , quickly 
evaporating and ejecting single and binary BHs. Thermal 
equilibrium with the surrounding stellar cluser is only re- 
stored when so few BHs remain that the subcluster's in- 
ternal timescale once again becomes commensurate with 
the cluster's interaction timescale. This process of segre- 
gation, decoupling, and evaporation has been extensively 
studied theoretically [SB], [26[ and numerically, both in full 
TV-body simulations (N « 10 3 - 10 5 ) 0, [H S EJ and 
in using approximations in larger N systems [la . Il8| . 

Previous studies have suggested that interactions in 
this dense cluster can produce many BH-BH mergers, 
whether through evaporated binaries [J] or throu gh run 
away BH-BH mergers in the dense cluster itself 3^, 31 
I32I ]. However, these studies have faced significant ob- 
jections to their details, because each has omitted some 
feature which could significantly inhibit the channel sug- 
gested. For example, the GW recoil produced during a 
BH-BH merger likely ejects the products of such a merger 
from any protocluster and prevents BH-BH merger run- 
away; the characteristic GW merger time of evaporated 
BH binaries depends sensitively (oc <t 8 ) on the inital clus- 
ter velocity dispersion; and even the initial cluster's dy- 
namics differs significantly depending on whether a dis- 
crete (two-component) or realistic (continuous) mass dis- 
tribution is used. 

For this reason, O'Leary et al. [HI performed a broad 
range of detailed numerical simulations of the dynamics 
of the segregated core and its coupling to the surrounding 
cluster, incorporating several critical physical features of 
dense cluster dynamics (a realistic black hole mass dis- 
tribution; a range of BH-BH merger kicks that includes 
recoil speeds larger than the most extreme cluster es- 
cape velocities, as suggested by the most recent numeri- 
cal simulations of unequal-mass binaries 33j, [3J] ; a range 



of cluster velocity dipersions a from 5 to 20kms _1 , en- 
compassing the range seen in present-day globular clus- 
ters; and both three-body and four-body interactions). 
These simulations form the basis of our estimates of the 
numbers and properties of the evaporated BH-BH binary 
population (A^m ~ 0.07A?bh) that is ejected during the 



decoupled phase, consistent with theoretical expectations 
25l 26j |. These "hot" ejected binaries have a thermal ec- 
centricity (e) distribution (2ede) and a lognormal binding 
energy (Eb = Gm\mil^a) distribution, with half the bi- 
naries having binding energies between 10 35 and 10 45 
times the thermal energy of the core, kT — (m*)cr 2 /3 
(see Fig. 6 of O'Leary et al.). Once ejected, the BH-BH 
binaries' orbits decay through GW emission, according 
to Peters' Eqs. (5.6-7) [35j], until eventually each binary 
merges. The initially most strongly bound and eccentric 
orbits merge first, followed by the initially widest and 
most circular [Fig. [Q, with a GW delay time t gw oc a 4 . 
Furthermore, because of the broad range of ejected bi- 
nary energies, while the differential distribution's char- 
acteristic merger time depends sensitively on the as- 
sumed cluster velocity dispersion (t gw oc a 4 oc c -8 ) and 
to a lesser extent on the precise details of the energy 
spectrum of ejected binaries, the delay time distribu- 
tion between 0.1 and 13 Gyr is comparatively robust, 
varying by ~ 50%. Specifically, when a population of 
rribh = 14M Q black holes evaporates from inside a clus- 
ter with a (T c i = lOkrns -1 velocity dispersion, the binary 
black holes will have a probability P(< t) of merging by 
time t, with 
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where we simplify the distribution of binding energies 
by assuming that it is flat in the log; a similar expres- 
sion applies to general velocity dispersions a c \, when 
the logarithm (only) is rescaled according to \ogt — > 
log£(cr c i/10 kms -1 ) 8 . Except for very early and very 
late times, this distribution is well approximated by 
dP/dt w 0.054/f, as expected from the relation between 
GW delay time (t gw oc a 4 ) and the logarithmic distribu- 
tion of semimajor axes implied by cluster evaporation. 
This simple 1/t decay in merger rate per cluster appears 
in the simulations of O'Leary et al. (see their Figure 5). 

Combining the number of black holes we expect per 
unit mass (Nbh ~ 3 x 10 _3 (M C //M Q )), the fraction of 
those black holes which should be ejected as binaries 
(A^bin ~ 0.07iVbh, based on a variety of simulations [l8(), 
and the rate of mergers per unit time given by Eq. |1} , 
we expect the number of cluster mergers per unit cluster 
mass after a time t to be TL c i{t) = 2.1 x 10~ 4 A/q 1 dP/dt. 
Assuming that clusters are formed from a fraction g c \ of 
available star formation and that only a fraction <7 CV ap of 
all cluster-forming mass possesses the birth conditions 
necessary for this process to occur, then the BH-BH 
merger rate R eV a.p(t) per unit comoving volume is 
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where dp/dt is the observed star formation rate per unit 
volume in the universe [361 . 37 1. While observations 
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FIG. 1: The distribution of GW merger timescales for circu- 
lar (solid step function) and eccentric (solid smooth function, 
based on 10 4 samples) 14Mq + 14M© binaries, assuming the 
binding energy is distributed logarithmically from 10 3 fcf,T to 
10 J fc;,T for fef,T the characteristic cluster energy. This figure 
assumes k b T = 0.5M Q (10kms- 1 ) 2 /3. 
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FIG. 2: The distribution of binary birth times for those binary 
black holes merging at present, given three different assump- 
tions for the star formation history in the universe: Nagamine 
et al. [37l ] (solid), Porciani and Madau's Eq. (5) (dotted), and 
(6) (dashed) 38]. In all three cases, clusters formed more than 
5 Gyr ago produce more BH-BH mergers than more recent 
cluster formation. 



have yet to converge to a unique star formation history, 
most models for dp/dt imply a present-day merger rate 
within 30% of « g ovap g c iMpc~ 3 Myr _1 ; they also univer- 
sally imply these binaries are born equally often in the 
nearby and early universe [Fig. |2J. Finally, based on 
the weak scaling of merger rate with a c \ implied by Pe- 
ters' equations, the BH-BH merger rate depends weakly 
(empirically, roughly logarithmically) on the assumed ini- 
tial cluster velocity dispersion <r c \ between 5km/s and 
20km/s. In short, binary BHs merging due to this process 
do so at a relatively well-defined rate, agreeing to order 
of magnitude with the merger rate expected if all mass 
presently in stars formed 10 Gyr ago: i? e vap ~ <7ci^eva.r,2 x 
IQ-H'I^pv/ilOGyr) 



Unless suppressed strongly (i.e., g c ig cvap <C 10~ 2 ), the 
BH-BH merger rate due to clusters will significantly ex- 
ceed the average rate densities expected from isolated 
stellar evolution, » 10- 2 Mpc~ 3 Myr _1 0, [39|. Addi- 
tionally, as noted in O'Leary et al., the BH-BH bina- 
ries produced from dynamical cluster evaporation will 
strongly favor pairs of the highest-mass binaries (e.g., 
14M Q + 14M or even 20M© + 20M Q ). On the other 
hand, isolated binaries rely on SN kicks and mass trans- 
fer to bring them to merging; therefore these merging 
BH-BH binaries are typically significantly less massive. 
Because merging BH-BH binaries produced in clusters 
possess a distinctive mass signature and could occur at 
unusually high rates, GW observatories can directly con- 
strain or even measure g c i3cvap- For example, in an op- 
timistic case (g c \ — g cvap = 1), the initial LIGO network 
could detect roughly 10 events per year, based on an 
estimated network range to 14M© + 14M Q binaries of 
125Mpc; the "enhanced LIGO" upgrade, with roughly 
twice the sensitivity, should see roughly 2 15 / 6 as many 
sources; and with roughly 20 x the range, the advanced 
LIGO network could detect as many as 3 x 10 4 (7 c i<7 e vap 
sources per year, permitting exquisite probes of early- 
universe cluster dynamics if indeed this process yields 
more BH-BH mergers than isolated stellar evolution. 55 1 

Gravitational wave observatories provide useful infor- 
mation precisely because g c \ and <? vap are so weakly 
constrained electromagnetically. For example, based on 
galactic @, Q and extragalactic @, H, EH cluster obser- 
vations, the fraction of stars born in clusters g c \ could 
be anywhere from 20% to 100%. Similarly, observations 
of galactic and extragalactic clusters cannot rule out all 
clusters more massive than 3 x 10 4 M Q and thus with more 
than 100 BHs undergoing runaway segregation and evap- 
oration. While such clusters are exceedingly rare in num- 
ber (p(M)dM oc Af- 2 HUSHED), they likely contain a 
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significant fraction of all cluster-forming mass: assuming 
that clusters range in size from 30M Q to lO 7 Af , roughly 
45% of all cluster-forming mass lies in these most massive 
clusters. While this fraction depends very weakly (loga- 
rithmically) on the limiting masses assumed for the clus- 
ter mass spectrum, a flatter mass function, as suggested 
by some observations 4^, 43 1 , would imply a significantly 
higher fraction of mass in these most massive clusters. 
Therefore, if all sufficiently massive clusters undergo seg- 
regation runaway, <? e vap could be close to unity. On 
the other hand, these evaporated BH-BH binaries appear 
only if their clusters of origin survived long enough as 
bound objects for gravitational segregation to occur. Be- 
cause gravitational segregation should occur much more 
rapidly than the age of observed long-lived clusters, such 
as the globular clusters of the Milky Way, g ova p must be 
larger than the corresponding fraction of a galaxy's mass: 
5evap ^ 10~ 4 — 10 -3 [lH- On the other hand, gravita- 
tional segregation should at best compete with and more 
likely occur more slowly than "infant mortality," the ten- 
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dency of roughly 70 — 90% of young clusters to disrupt 
within their first ~ 10 Myr due to photoionization- and 
SN -driven gas ejection [E EE EE EE EE E3 ■ Based on 
only 10% of all clusters surviving a rapid "infant mortal- 
ity" epoch and on 45% of all clusters being sufficiently 
massive for the Spitzer instability to occur, we expect 
5ovap ~ 5 x 10~ 2 . To summarize, we expect g c i<?evap 
could be as high as 1 (corresponding to 10 merger de- 
tections per year with initial LIGO); is likely 5 x 10~ 2 
(1 event every two years with initial LIGO; 3 events per 
year with "enhanced" LIGO); and is very likely higher 
than 10 -4 (producing merging BH-BH binaries slightly 
less frequently than isolated binary stars, but with sys- 
tematically higher masses). 
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